Mechanical relaxation in glasses and at the glass transition 



U. Buchenau 

Institut fiir Festkorperforschung, Forschungszentrum Jiilich, Postfach 1913, 52425 Jiilich, Germany 

(February 1, 2008) 

The Gilroy-Phillips model of relaxational jumps in asymmetric double well potentials, developed 
for the Arrhenius-type secondary relaxations of the glass phase, is extended to a formal description of 
the breakdown of the shear modulus at the glass transition, the a-process. The extension requires 
the introduction of two separate parts of the barrier distribution function f{V), with a different 
temperature behavior of primary and secondary parts, respectively. The time-temperature scaling 
of the a-process, together with a sum rule for the whole barrier distribution function, implies a strong 
rise of the integrated secondary relaxation with increasing temperature above the glass transition. 
Thus one gets for the first time a quantitative relation between the fragility of the glass former and 
the fast rise of the picosecond process observed in neutron and Raman scattering. 

The formalism is applied to literature data of polystyrene, vitreous silica and a sodium silicate 
glass. In the glass phase of polystyrene, one finds a temperature-independent secondary barrier 
distribution function, in agreement with an earlier Raman result from the literature. Above the glass 
transition, the secondary barrier distribution function increases with temperature as predicted. The 
findings allow for a new interpretation of the fragility and the entropy crisis at the glass transition. 



I. INTRODUCTION 

Relaxation in glasses, sometimes also called secondary 
relaxation to distinguish it from the primarv. relaxation 
at the glass transition, is generally believedETtl to be well 
described in terms of the Arrhenius-Kramers pictured, 
with a relaxation time ry given by the Arrhenius rela- 
tion 



TV = Toe^/^- 



(1) 

where rg is a microscopic time of the order of 10^ ^'^ sec- 
onds, V is the energy of the barrier between two energy 
minima of the system, and T is the temperature. 

In contrast, the primary relaxation or a-process, the 
onset of the flow process at the glass transition tempeE|u. 
ture Tg and above, seems to follow a much steeper lawQ'El 



Toe 



A/(T-To) 



(2) 



where A and Tq are constants with the dimension 
of a temperature. This is the well-known empirical 
Vogel-Fulcher-Tamman (VFT) or Williams-Landel-Ferry 
(WLF) equation. Tq, the Vogel-Fulcher temperature, is 
smaller than Tg] the closer it lies to Tg, the more fragile 
is the glass former. 

Since the, Arrhenius law has a sound microscopic 
backgroundD and the VFT or WLF equation has not, it 
seems reasonable to build a joint quantitative description 
on the former, bearing in mind the physical difference of 
the two processes. This is the intention of the present 
paper. 

In glasses, one has to reckon with a whole distribu- 
tion of relaxational jumps, not only over different poten- 
tial barrier heights, but also between energy minima of 
different energy. Thus one has to extend the classical 



Arrhenius-Kramers treatment! of a thermally activated 
relaxation process in a symmetric double-well potential 
to deal with a broad distribution of barrier heights and 
asymmetries. A distribution in the barrier heights was 
first considered by FrohlichQ, but only for symmetric po- 
tentials. An asymmetric multiminimum situation was 
solved by Hoffman and Pfeiffem. But none of these early 
attempts dealt simultaneously with a broad distribution 
in both quantities, the barrier height and the asymmetry 
of the wells. The necessity of such a double distribu- 
tion was first recognized inJ;he tunneling model for the 
two-level states bektw 1 Kla in 1971. Ten years later, 
Gilroy and PhillipsE£l extended the scheme to classical 
relaxation processes. Also, they were the first to draw 
a parallel between mechanical relaxation data and the 
quasielastic part of the Raman scattering. At the end of 
the same decade, the soft potential model postulated a 
relation between the tunneling states arid, the low barrier 
classical relaxation processes in glassesEJ. Nevertheless, 
up to now onlyjiejiL checks of theae-nr^tulates for cbojiami- 
cal mechanicalBEra, RamanElOt3 and neutronB data 
have been reported in the literature. 

The present paper begins in Section II with a discus- 
sion of the Gilroy-Phillips model, and a derivation of its 
connection to rheology. It turns out that one can define 
a barrier distribution function f{V) to describe the me- 
chanical shear relaxation at different temperatures and 
frequencies. The integral of this barrier distribution func- 
tion over all barrier heights V must equal 1 to bring the 
shear modulus down to zero. In order to include the flow 
process into the same scheme, one separates f{V) into 
two parts, fs{V) and fa{V). The first of these describes 
the secondary relaxations in the glass phase, the second 
describes the a-process in the undercooled liquid, respec- 
tively. As will be seen, the sum rule for the total barrier 
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distribution function supplies a quantitative basis for An- 
gell's conjectureEj of a relation between the fragility and 
the rise of the fast picosecond process above the glass 
transition. The low-barrier part of fs{V) determines the 
relaxational part of neutron and Raman scattering. The 
corresponding equations are derived. 

Section III apphes the equations to determine f{V) 
from literature data for amorphous polystyrene, vitre- 
ous silica and a sodium silicate glass. From the results, 
one gets a first impression as to whether one gets the 
same secondary barrier distribution function from differ- 
ent methods, in particular if one compares the high fre- 
quency neutron, Raman and Brillouin scattering results 
with the low frequency torsion pendulum or creep data. 
Furthermore, one gets a feeling for the amount of reduc- 
tion of the shear modulus by the secondary relaxation 
processes in the glass phase. As we will see, the findings 
suggest a generalized Maxwell criterion for the onset of 
the glass transition, namely that the flow begins when 
the shear response from the secondary relaxation equals 
the elastic one. Section IV compiles and discusses these 
results and their possible significance for our view of the 
glass transition. Summary and conclusions are given in 
Section V. 



II. THE GILROY-PHILLIPS MODEL 



The interaction between the shear strain and the De- 
bye relaxation center is described by the change of the 
asymmetry A under the influence of the strain. The in- 
teraction is characterized by the coupling parameter 7, 
leading to an asymmetry A -I- 76 of the relaxation in the 
strained glass. 7 must be considered to depend both on 
V and A. 

The free energy F of the relaxation center reads 



F = -kgT hi 



2 cosh 



A- 



-76 



2kBT 



(3) 



which has the second derivative with respect to the shear 
distortion e 



d^F 



7 



4fcBTcosh^(A/2fcBT) 



(4) 



The second derivative determines the contribution of that 
specific relaxing entity to the difference between the shear 
moduli at infinite and zero frequency. The equation 
shows that the main influence on the shear modulus is 
due to relaxation in potentials with asymmetries smaller 
than ksT] for larger asymmetries the influence decreases 
rapidly because of the square of the hyperbolic cosine in 
the denominator. 



A. The asymmetric double- well potential 



B. The barrier distribution function fiV) 



Let us denote the shear strain by e, the shear stress by 
<T and the (infinite frequency) shear modulus by G. G 
will generally depend on the temperature T. 




FIG. 1. Asymmetric double well potential with barrier 
height V and asymmetry A as a function of a generalized 
coordinate. 

The structural relaxation is taken to be a superposi- 
tion of independent Debye relaxation centers in asym- 
metric double- well potentials with two minima, as shown 
in Fig. 1. The energy of the left minimum is —A/2 and 
the energy of the right minimum is -f A/2. The height of 
the barrier is V . 



We want to calculate the frequency dependence of the 
shear modulus under the assumption of slowly varying 
distribution functions in the parameters V and A. In 
detail, we assume a number density of relaxing entities 
n(y. A) and a coupling constant 7(y, A) which are both 
approximately constant if either or A is varied by an 
amount of the order of the thermal energy ksT. 

Under this assumption, it is safe to neglect as well the 
influence of the asymmetry on the relaxation time. We 
assume the relaxation time tv to be given by the Arrhe- 
nius equation 

We then integrate over the asymmetry A to obtain the 
step 5G between the shear moduli at infinite and zero 
frequency from all relaxation centers with barrier height 
between V and V + dV 



SG = dV 



^ 7^n(F, A)dA 

ex. 4A;BTcosh^(A/2fcBT) 



(5) 



Since one has only contributions in the near neighbor- 
hood of A = 0, where n{V, A) « n{V, 0), and since 



dA 



cosh^(A/2fcBr) 



one finds 



dG = "f^n(y, 0)dV. 



(6) 



(7) 
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This is different from a single relaxation in a symmet- 
ric potential, where the step in the modulus is inversely 
proportional to the temperature. The physical reason for 
this difference is clear: As the temperature rises, relax- 
ation centers with higher and higher asymmetry begin 
to contribute to the step in the modulus. This is an 
important difference between relaxation in crystals and 
relaxation in disordered matter. 

The temperature-independent step in the modulus is 
determined by the barrier distribution function f{V), de- 
fined by 



f{V) 



G 



(8) 



This parameter combination can be argued to remain in- 
dependent of temperature, even if G varies with temper- 
ature, considering the relaxiftg entity as a small misfit 
region in an elastic mediumEd, a misfit region which is 
able to change the sign of the misfit by jumping over the 
barrier. Here, however, this argument will not be given 
in detail. 

The frequency dependence of the complex shear mod- 
ulus at the frequency u and the temperature T reads 



G'{u;,T) _ G 
G ^ G 



= fiv) 



dV 



G 



fiV) 



UTydV 

1 + uj^tI 



(9) 



(10) 



where Ty is a function of V by the Arrhenius relation eq. 
(|l|), and Ge is the zero frequency modulus after the decay 
of all the relaxations in the system. 

These two equations describe the real and the imagi- 
nary part of the frequency-dependent shear modulus at 
all frequencies and temperatures. As long as one can 
reckon with a temperature-independent number of un- 
coupled relaxation centers, the barrier distribution func- 
tion f{V) remains temperature- independent. 

But even if f{V) depends on temperature, it can still 
be formally used to characterize the shear relaxation. As 
shown in the next subsection, cp£ can rewrite the con- 
ventional rheological expressionscj in terms of f{V). The 
advantage of the choice of f{V) lies in the possibility 
to distinguish the trivial Arrhenius temperature depen- 
dence from other, nontrivial temperature changes. These 
nontrivial temperature changes will then reflect in a tem- 
perature dependence of f{V). 



C. Rheological equations in terms of f{V) 

Comparing the two expressions, equs. ^ and (0), to 
those in the textbooks, for instance the one on polymers 
by Ferrycj (chapter 3, equs. (23) and (24)), one finds 
the relation between the rheological relaxation function 
H{t) and the barrier distribution function f{V) 



if(Toe^/'=-^) = H{tv) - GkBTfiV). 



(11) 



With this equation, one can rewjrite all the exact and 
approximate rheological relation£3 in terms of f{V). To 
do this, one first has to define a convenient equivalent 
function 1{V) to the conventional rheological function 
L{t), which is needed whenever one wants to calculate a 
compliance 



L(Toe 



V/keT 



) = Liry) 



kBTljV) 
G ' 



(12) 



In the following, thet4iiost important equations of chap- 
ter 3 of Ferry's bookcS are translated into the Gilroy- 
Phillips notation. Ferry's equation (19) for the time- 
dependent modulus G{t) reads 



/oc 
iJ(T)e-*/^dlnT 
-OO 

and translates into 

roc 

Git) ^Ge + G / /(F)e-*/^^dF. 



(13) 



(14) 



For a viscoelastic liquid, the zero frequency modulus 
Ge = 0, so one must have 



fiV)dV = 1. 



(15) 



This is the sum rule for the barrier distribution function 
f(V). It has important consequences for the connection 
between primary and secondary relaxation, as discussed 
in the next subsection. 

Next, there is Ferry's eq. (20) for the compliance 

1 t 

J{t) = -+ L(r)(l-e-*/^)dlnr+-, (16) 
G J-oc Vo 

where rjo is the viscosity, which translates into 



Jit) 



1 

G 



l{V)il 



-t/rv 



)dV 
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Vo' 



(17) 



The viscosity r/o can be calculated from Ferry's eq. (28) 

/OO 
TH{T)d\nT, (18) 
-OO 

which translates into 

poo 

Vo = G Tvf{V)dV. (19) 
Jo 

Ferry's equs. (21) and (22), the transformation from 
H{t) to L{t) and back, read 



H 



t/u — 1 



dim 



(20) 
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and 



H = 



L 



J. + f"' 

G J—oo 1—u/t rio 



— U/T 

They translate into 
1{V) = 

with 

IfiV) 

and 

f{V) = 



(21) 



f{E)dE 
exp((F - E)/kBT) - 1 



1 



TvG 



h{V) +{TTkBTl{V)Y 



(22) 



(23) 



(24) 



with 



Ii{V) 



l{E)dE 
l-exp((^-y)/fcBT)- 



(25) 



With these exact equations, one can calculate the me- 
chanical response for any type of shear experiment for a 
given barrier distribution function f{V). The reverse, the 
determination of fiV) from experimental data, is more 
difficult, because the exact equations are integral equa- 
tions. Nevertheless, one can start to determine a first 
approximation to f{V) from measurements of G' and G" 
using the crude approximations 



G 



f{V)dV 



and 



/(fcsrin(lMo)) ^ - - ^" 



nGkeT 



(26) 



(27) 



D. Primary and secondary relaxation 

It is quite clear that one needs to distinguish secondary 
and primary processes, because their physical mechanism 
is different. Thus one has to distinguish between fs{V), 
the secondary barrier distribution function of the sec- 
ondary Arrhenius relaxation, and fa{V), the primary 
barrier distribution function for the primary a-process 
or flow process (see Fig.r^). 

As is well established^^, Hair/Ta) is independent of 
the temperature. This is the time-temperature scaling 
of the a-process, sometimes also denoted as thermorhe- 
ological simplicity. Since the factor GksT in eq. (|Tl] ) 
between H{t) and f{V) varies only weakly with temper- 
ature, this implies that the primary barrier distribution 



ftinction ( V) is an essentially temperature-independent 
function of — Va{T), where Va{T) denotes the maxi- 
mum of this strongly peaked function. From the Vogel- 
Fulcher law eq. (^, one expects the temperature depen- 
dence 



V^{T)^V^iT,) 



T{Tg ~ To) 



(28) 



showing the divergence of the fictive Arrhenius barrier of 
the flow process towards the Vogel-Fulcher temperature 
To. 
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FIG. 2. Secondary barrier distribution function fs{V) (the 
shaded area), together with the cutoff by the a- relaxation 
peak (schematic). 

If time-temperature scaling holds, the weight Wa of the 
a-process, given by 



Wo 



UV)dV 



(29) 



should be temperature- independent. In the comparison 
to experiment, we will see that Wa tends to be close to 
1/2. 

What does this imply for the secondary relaxations? 
The a-process is also an upper cutoff for the secondary 
relaxation; at the end of the process, the long time shear 
modulus is zero. If a secondary relaxation barrier is too 
high, the relaxing entity will flow away before it has a 
chance to jump. Therefore there is a steep cutoff for the 
secondary barrier distribution function fsiV) at Va(T). 
According to the sum rule eq. (Esh 



K,(T) 



fs{V)dV = l-wa^ const « i. (30) 



With this sum rule, the decrease of Va{T) with in- 
creasing temperature implies that the secondary barrier 
distribution function fsiV) must increase with increasing 
temperature. This increase will be stronger for more frag- 
ile glass formers. Such a connection between the fragility 
and the rise of the fast process above Tg has been indeed 
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postulated empirically@; here we will be able to quantify 
this connection. 

The increase of fs{V) above the glass temperature can 
be characterized to first order by 



f,{V) = fsiV,Tg) 



l + asiV) 



(31) 



by 



One can define an average temperature coefficient as 



1-Wa 



as{V)fs{V)dV. 



(32) 



Differentiating the sum rule eq. ( |30| ) for fs[V) with 
respect to the temperature, calculating the derivative of 
Va (T) with respect to temperature from eq. ( p8| ) and us- 
ing eq. (|3^), one finds the relation between the fragility 
and the average rise of the secondary relaxation 



(33) 



This relation will become clearer in the next subsec- 
tion, where the simplest possible case of a constant sec- 
ondary barrier distribution function is discussed. 



E. The generic case fs(V) = const 

The deep implications of the Gilroy-Phillips formula- 
tion of the a-process are more clearly seen in the sim- 
plest possible case, shown in Fig. 2. Let us assume 
fs{V) = const and as{V) — const. If the rise of fsiV) 
with temperature is strictly linear, equs. ( ^o|) and (|3ll ) 
imply 



Vo.{T) 



Tg + asiT-TgY 



(34) 



which is not exactly equal, but very close to the empirical 
Vogel-Fulcher-Tamman or Williams-Landel-Ferry equa- 
tion (S) , with the Vogel-Fulcher temperature Tq given by 



Tn = r„ 1 - 



1 



(35) 



It is obvious how this comes about: at Tq, the den- 
sity of secondary processes extrapolates to zero. Thus 
one has to proceed to infinitely high barriers to satisfy 
the sum rule (^0|). This gives a completely new view 
on the puzzling fragility of glass formers: the abnormal 
temperature dependence is a consequence of the time- 
temperature scaling of the a-process, and of a strictly 
linear rise of the number of secondary relaxing units with 
temperature. 

The decrease of the secondary barrier distribution 
function fsiV) implies a decrease of the number oiopin- 
ima of the glass former. Thus one gets an equalityn3 be- 
tween the Vogel-Fulcher and the Kauzmann temperature 



Tk, the latter being defined as the temperature where the 
excess entropy of the glass former over the corresponding 
crystalline system extrapolates to zero. If there are no 
minima between which the glass former can jump, there 
is no excess entropy, the old Adam-Gibbs ideaEJ. 

Note that eq. (^5|) for the Vogel-Fulcher or Kauzmann 
temperature holds not only in the generic case of a con- 
stant secondary barrier distribution function, but for any 
fs(y), as long as one can reckon with the same strictly 
linear temperature rise of the function for all V. 

To complete the discussion of the glass transition 
peculiaritiesQ in the Gilroy-Phillips picture, let us look 
at the stretching of the a-process, empirically described 
by the Kohlrausch equation 



G(t) =G„ exp(-(T/T„)'^ 



(36) 



where the Kohlrausch exponent f3 liesH between 0.3 and 
0.7, and Gq is a free parameter. The smaller f3 is, the 
more stretched is the a-relaxation, and the stronger it 
deviates from a simple exponential decay. 

In order to calculate /3, the definition of the primary 
barrier distribution function faiV) must be more spe- 
cific. Let us assume a gaussian centered at Va with 
weight Wa — 1/2 and with a full width at half maxi- 
mum (FWHM) denoted by Aq,. Calculating G{t) from 
eq. (ttq), one then finds fairly Kohlrausch- like curves for 
the simple generic case of Fig. 2, at least in the time 
region of the a-process, with Ga ~ 0.55. .0.75 G and Tq 
about 80 % of the Arrhenius value for Va ■ There is a de- 
viation of the Kohlrausch fits from the calculated curves, 
but it is so small that it would be hard to see in an ex- 
periment. It turns out that the ratio Aa/Va determines 
the Kohlrausch exponent /3; if it is 0.05, then /3 « 0.7; 
for the ratio 0.1, /? « 0.5 and for the ratio 0.2, f3 w 0.3. 
So the broader fa{V), the more stretched the relaxation, 
not unexpected. 

Note this is merely a change of description. The 
Gilroy-Phillips formulation does not really explain the 
puzzling features of the glass transition, the fragility, the 
entropy crisis and the stretching. But it supplies a de- 
scription which allows to look for a new kind of explana- 
tion. We will return to this point in the discussion. 

In order to measure a^iV) in the picosecond range by 
scattering methods, one still needs the equations for the 
neutron and Raman scattering functions in terms of the 
barrier distribution function. These will be derived in 
the next subsection, the last part of the description of 
the Gilroy-Phillips model. 



F. Neutron and Raman scattering 

One can carry out the same integrations over asym- 
metries and barrier heights as in the shear relaxation for 
the neutron scattering cross section. Let us begin with a 
single asymmetric double well, let us assume that atom 
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j has a coherent scattering length bj and an incoherent 
scattering cross section aj, and that it jumps from the 
position —dj/2 to dj /2, with the origin of the coordinate 
system in the middle between the two. 

For the incoherent jorielastic scattering in the one- 
phonon approximationEj, it suffices to calculate the mean 
square displacements. These can be obtained from the 
Boltzmann occupation factors of the two minima of the 
potential. The average position vector fj of atom j is 
given by 



d. , A 
— tanh — — — 
2 2kBT 



and its average square is 



^2 7 

< r >— — 



(37) 



(38) 



so the mean square displacement contribution of the re- 
laxation to atom i reads 



< r. 



>2= 



4 



4cosh^ A/2fcBT 



(39) 



In this expression, we recognize the same inverse cosh as 
in eq. (^), which can be again integrated over the asym- 
metries, if the jump eigenvector stays essentially the same 
for the different local asymmetries. 

The scattering contribution is a Lorentzian with a half 
width at half maximum in w given by the inverse relax- 
ation time, determined by the barrier height according 
to the Arrhenius equation (0). The weight of the contri- 
bution in S{Q, Lo) is determined by the number of relax- 
ations and the weighted sum of the jump vectors 



orientational averaging over the structure factor of the re- 
laxation. Even after this averaging, the structure factor 
need not show a simple behaviour like the incoher- 
ent one^but does still contain information on the jump 
vectorstj. 

The scattering measurements do not give f{V), but 
rather the product n(y^,0)d^. If one wants the propor- 
tionality factor between those two quantities, one needs 
additional information about the relaxing entities. How- 
ever, there is an elegant and general way to obtain this 
proportionality factor for very low barriers from the soft 
potential modelEj, which describes the tunneling states 
and the low-barrier classical relaxation as similar modes 
with a double-well potential distribution. It is not very 
difficult to derive an equation for the barrier distribution 
function f{V) as defined here in terms of the definitions 
in this paper. One finds 



2C 



^3/4yl/4' 



(43) 



where C (in principle Ci for longitudinal waves and Ct 
for transverse waves) is a dimensionless constant of the 
order of 10~*, which can be taken from acoustic attenu- 
ation measurements below 4 K, and W is the crossover 
energy between tunneling and vibrational modes, which 
can be measured from the crossover regions of the specific 
heat or the thermal conductivity at low temperatureal^. 
The soft potential model has a fourth parameter Pg for 
the density of these modes. With this parameter, the 
relation between n{V, Q)d^ and the secondary barrier dis- 
tribution function /s(T^) reads 



(44) 



(40) 



The sum is over all atoms in the sample, and a is their 
average incoherent cross section. 

Integrating over the barrier heights as well, one gets 
the equation for the incoherent scattering 



Sinc{Q,^^) = Vail 



fesTln— ,0 

LOTq 



klT^Q^d^ 
3w 



(41) 



where Va is the atomic volume, and the prefactor 1/3 
stems from the directional average. This equation is 
again an approximation, which holds if f{V) does not 
vary strongly with V. The coherent scattering is ob- 
tained replacing d"^ by d^^y^{Q) with 



dlohiQ) 



<|^6,e''3«^Qd;- |2>, 
j 



(42) 



where p is the mass density. For glasses consisting of 
more than one sort of atoms, the relation might fail if 
strongly and weakly scattering atoms have different jump 
widths. But with this relation, one can determine the sec- 
ondary barrier distribution function fs{V) from neutron 
scattering measurements without adaptable parawieter, 
using soft potential parameters from the literaturalj. 

The Raman scattering from relaxations in glasses is 
not so easily calculated. However, experienceEjE3 shows 
that neutron and Raman scattering give the same spectra 
as long as one stays at frequencies well below the boson 
peak. Thus one can use the Raman spectra as one uses 
incoherent neutron scattering data, with the disadvan- 
tage of an additional general adaptable parameter for the 
overall intensity. The advantage of the Raman technique 
is a much higher intensity and a much better resolution, 
allowing one to assess much lower frequencies. 



where Rj is the equilibrium position of atom j, 6 is the 
average scattering length and the brackets indicate the 



III. COMPARISON TO EXPERIMENT 
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A. Values from different techniques 



Let us first consider which barriers one samples with 
a given technique. Creep measurements cover the time 
range from 0.1 seconds to several weeks. This is a mea- 
surement in the time domain, applied mostly to mea- 
surements of the a-process at the glass transition, where 
one should not use the Arrhenius relation eq. (|l|). Nev- 
ertheless, one can formally calculate barrier heights of 
the order of 27 to 41 ksTg, around 1 eV for polystyrene 
(Tg — 373 K) and around 4 eV for vitreous silica 
{Tg = 1473 K). 

The torsion pendulum method with frequencies around 
1 Hz sees relaxations around 0.1 seconds. In terms of the 
Arrhenius relation with tq = 10"'^'^ s, this implies a bar- 
rier height of 50 meV at 20 K. In order to get to 1 eV, one 
has to rise the temperature by a factor of twenty, thus 
getting close to the glass temperature of polystyrene. 

Proceeding to higher frequencies, on has the vibrat- 
ing reed measurements around 10 kHz, ultrasonic data 
in the MHz range, light scattering Brillouin data around 
10 GHz and, finally, Raman and neutron data between 
a few GHz and a few hundred GHz. For the latter two, 
the lower limit holds only for the Raman technique; if 
one looks for the weak quasielastic scattering from sec- 
ondary relaxations, the neutron technique in practice has 
a lower limit of 100 GHz. The upper limit of about 300 
GHz is given by .thr crossover from relaxational to vibra- 
tional scatterin^lJo. Thus one sees only the uppermost 
frequency band of the relaxational scattering with neu- 
trons. Nevertheless, neutrons play an important role, be- 
cause they serve to validate the Raman scattering data. 

For the fast relaxation at 200 GHz, the relaxation time 
is of the order of a picosecond, only ten times longer than 
the microscopic time scale of the vibrational motion. The 
Arrhenius relation translates this again into a barrier of 
50 meV at 300 K. Thus neutron and Raman measure- 
ments at room temperature sample the same relaxations 
that one sees in a torsion pendulum measurement around 
20 K. If one is still in the glass phase at room tempera- 
ture, one can thus check the temperature independence 
of the barrier distribution function f{V) (in that case, 
f{V) = fs{V)) by a comparison between a torsion pen- 
dulum and a neutron or Raman experiment. Naturally, 
the same can be done by a comparison of a torsion pen- 
dulum and a Brillouin scattering experiment. 

If one wants to determine the barrier distribution func- 
tion from the mechanical damping for different frequen- 
cies, one often has to compare measurements of differ- 
ent elastic constants. A torsion pendulum measurement 
provides immediately the real and imaginary part of the 
shear modulus G. The vibrating reed technique measures 
Youngs modulus Y. In terms of the bulk modulus B, the 
inverse of the compressibility. Youngs modulus reads 



Since that quotient is close to 3G, it is a reasonable ap- 
proximation to identify the tan (5 = G' /G" results of the 
vibrating reed technique (in that technique often denoted 
as = tan (5) with that of the shear modulus. The real 
part of the shear modulus can then be estimated to be 
Y'/3, or calculated accurately from the above equation 
if the bulk modulus is known. 

A less safe connection is to longitudinal sound mea- 
surements, which determine the real and imaginary part 
of the elastic constant cn, sometimes also denoted as the 
modulus M. In terms of B and G, cn reads 



Cll 



4 

B + -G. 



(46) 



Both contributions are of comparable size. Therefore, if 
one wants to determine tan 6 for the shear modulus from 
longitudinal sound attenuation data, one has to know 
the relative size of the two damping contributions from 
compression and from shear. Here it should be noted 
that the mechanical damping at very low temperatures 
is quite similar ieir transverse and longitudinal phonons 
in many glassesES. If one believes the soft potential hy- 
pothesis for a common origin of tunneling states and low 
barrier relaxation in glasses, the same should hold in the 
temperature range up to 50 K. Then one can take tan (5 
from the inverse mean free path l^^ of the longitudinal 
sound wave according to 



tan 5 ~ 



vl- 



2uj ' 



(47) 



where v is the sound velocity. 

It is better to use data from transverse sound waves; 
there, the damping relates directly to the shear. How- 
ever, in particular for light scattering Brillouin experi- 
ments, it is much easier to determine the damping for the 
longitudinal waves, which provide a much clearer signal. 
That damping is often given in terms of the half width 
at half maximum F of the Brillouin line. Then 



F 

tano = — , 
2u;' 



(48) 



Y 



9)BG 
3B + G' 



(45) 



where lo is the frequency of the Brillouin line. 

Finally, one problem of the determination of f{V) is to 
know G at the temperature of the measurement; usually 
one has only G" and G" . A way out of this problem is 
to start at low temperatures, where G w G', determine 
f{V) via eq. (p7|), and then pursue G'/G to higher tem- 
peratures via eq. (|26|). This way was followed throughout 
in this paper, at least as far as the determination of the 
secondary relaxation was concerned (the primary relax- 
ation is too sharp to use the crude approximation of eq. 
(p7|); there, one has to use the exact equations). 

One can try to check the G(T) values by Brillouin light 
scattering measurements of the transverse sound waves, 
usually at frequencies around 10 GHz. Even there, one 
has to reckon with some influence from the low-barrier 
part of the relaxations; this is, however, only a small cor- 
rection which can be easily done if one has f{V), using 
eq. (H). 
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B. Polystyrene: secondary relaxation 

We begin the comparison to experiment with a heav- 
ily studied glass formet, atactic polystyrene, one of the 
most fragile substanceaj, an amorphous polymer where 
one can rely on a large number of experimental data, 
both at low temperatures and at the glass transition. 

The low temperature data were evaluated for temper- 
atures above 10 K; at that limiting temperature, one 
can begin to reckon with the validity of the Kramers 
picturey. Fig. 3 shows a compilation of many mechan- 
icalJow tempepature data: torsion pendulum data at 1 
H^ZL at. 6 HzHa, vibrating reed data at 3, 34 and 87 
kHzEjtil and BtilJouin damping of longitudinal sound 
waves at 10 GHa^. The data cannot be said to coincide 
perfectly in this Gilroy-Phillips evaluation; nevertheless, 
the agreement is fair enough to support the concept of 
a constant number density of uncoupled relaxation cen- 
ters. There is no systematic variation with frequency; 
one rather has the impression that the differences stem 
from the different sample preparation of these seven mea- 
surements. The line in Fig. 3 represents the fit of fs{V) 
to these mechanical data. It is seen that fs{V) rises to- 
wards low barriers, as one would expect from the soft 
potential equation (p3[). But the fit line falls below the 
soft potential expectation already at rather low barriers, 
similar to observations in other glassestil. 
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FIG. 3. Secondary barrier distribution function fs{V) cal- 
culated from literature data of the mechanical damping of 
amorphous polystyrene in the glass phase at different frequen- 
cies. References see text. The line is a fit; the same fit is also 
shown in Figs. 4, 5 and 6. 

Fig. |-Ut compares the same fit to the evaluation of 
RamanO and neutronEj data in terms of eq. (^). In 
the Raman case, data were adapted by an appropriate 
multiplication factor. The neutron data were evaluated 
using eq. (^3) and the soft potential fit parameters of 
polystyreneEiy without any free parameter. The fit of 
the neutron data requires a good description of the vi- 
brational part of the scattering; details will be given in 



a forthcoming publication combining time-of-fiight and 
backscattering data of polystyreneEj. 

The good agreement between data points and the me- 
chanical data fit line in Fig. 4 corroborates the earlier 
conclusiorJlZI of a temperature- independent fs{V) in the 
glass phase of polystyrene. Note that the earlier con- 
clusion does not stem from a comparison of Raman and 
mechanical data, but rather from a comparison of Ra- 
man data at three different temperatures, namely 100, 
200 and 300 K. The two different ways to check the tem- 
perature behaviour of fs{V) in the glass phase provide 
the same result. 
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FIG. 4. Secondary barrier distribution function fs{V) cal- 
culated from literature data of neutron and Raman scattering 
from amorphous polystyrene in the glass phase. References 
see text. 
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FIG. 5. Secondary barrier distribution function fs{V) cal- 
culated from literature data of neutron and Brillouin data of 
amorphous polystyrene, both in the glass phase and above 
the glass temperature Tg. References see text. 

This temperature-independence, however, does no 
longer hold in the undercooled liquid phase, above the 
glass temperature of 372 K of polystyrene. Fig. 5 shows 
neutrorS and longitudinal smmd wave damping data 
from the Brillouin techniqueo above Tg. In order to 
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relate to the preceding figures, f{V) is again plotted 
against the barrier height V. This implies that one sees 
the onset of the increase of fs{V) with increasing tem- 
perature at different values of V in the two techniques, 
at 0.074 eV for the neutrons and at 0.163 eV for the 
Brillouin data. Note that in both cases the frequency is 
too high to see the a-process at the temperatures of the 
measurement. 

We conclude that fs{V) does indeed increase above 
Tg, as postulated above on the basis of the sum rule for 



f(y), eq. (|15|), and on the basis of the temperature de- 
pendence of the a-process. The rise of fs{V) above Tg 
can be characterized by the linear relation eq. (^ij) with 
as — 7.7 ± 1 for the neutron data and as — 10 ± 3 for the 
Brillouin data (in the latter case, the large error is due to 
the small number of points and the insecurity of the value 
at Tg), within experimental error the same temperature 
coefficient for both sets of data. 

Fig. 6 shows fsiy) for polystyrene over the whole bar- 
rier range, together with the results oLa torsion pen- 
dulum measurement at 1 Ha up to TgEZI, an ultrasonic 
measurement just below TgL3 and the neutron result up 
to TgEj. The shaded area represents fsiV); the a-peak 
above 1 eV will be discussed in the next subsection. 




0.4 0.6 0.8 
barrier height V (eV) 

FIG. 6. Secondary barrier distribution function fs{V) cal- 
culated from literature mechanical damping and neutron data 
of amorphous polystyrene up to the glass transition. The peak 
at the end shows the gaussian fa{V) describing the a-process 
at the glass temperature Tg. References see text. 



In the case of polystyrene, one finds Wa = 0.54 for the fit 
function in Fig. 6. 

The width — 0.107 eV, the glass temperature 
Tg = 372 K and iJafi coefficient as = 6.9 ±0.1 was 
fitted to creep dataO of polystyrene with a molecular 
weight of 600000 g/mol, data which are also treated in 
Ferry's bookE3 and which are shown in Fig. 7. The 
fit requires as an additional parameter the infinite fre- 
quency modulus G = 1.69 GPa at 373.8 K. Together with 
the low temperature data and the known density varia- 
tion of polystyrene with temperaturecH, one calculates 
a Griineisen Tg — 3.8 for G in polystyrene. With this 
value, one can also calculate G for higher temperatures. 
The creep function J(t) is obtained by first calculating 
1{V) via eq. ( ^2|), the viscosity from eq. ( p^ and finally 
J(t) from eq. (\17\). 




FIG. 7. Creep data of the a-process in polystyrene at the 
glass transition (reference see text), together with a fit in 
terms of f{V). 

As seen in Fig. 7, one can describe the temperature 
shift of the a-process with an appropriate rise of fs{V) 
above Tg. The sum rule eq. ( pO|) forces a temperature 
shift of Va (T) , which in turn provides the experimentally 
observed shift factors towards higher temperatures. 

The scheme works quite well up to three decades in 
compliance; for still higher compliances, one sets into the 
plateau regime from the chain cntangiementtS, which is 
beyond the present considerations. 



C. Polystyrene: glass transition 

The a-process is characterized by the position Va (T) , 
the full width at half maximum Aq, of faiV) and the 
weight Wa. The latter two should be temperature- 
independent. Let us define Va{Tg) — 15fcBTgln(10) in 
order to have the corresponding relaxation time at 100 
s, and let us choose a gaussian for the primary barrier 
distribution function fa{V). 

If one knows fs{V) from measurements in the glas s 
phase, Wa can be calculated from the sum rule eq. (pO|). 



D. Vitreous silica and sodium silicate 

Vitreous silica is the case for which the Gilroy-Phillips 
ideaE3 was originally developed. The compatibility be- 
tween mechanical relaxation around 20 K and inelas- 
tic neutron scattering at roam temperature was demon- 
strated seven years lateiO. A very recent Raman 
experimentta showed the temperature independence of 
fs{V) up to room temperature, together with an excellent 
agreement of the shape of the function from mechanical 
data with the one from Raman data. 
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FIG. 8. The barrier distribution function of vitreous silica, 
together with torsion pendulum and Brillouin damping data. 
References see text. 
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FIG. 9. Dynamical mechanical data of vitreous silica at the 
glass transition, together with the model fit (see text). 

Of the many low-barrier mechanical data in the litera- 
ture, we include here merely one of the Brillouin damping 
experiments at 35 GHzlj. The range of higherj-barriers 
is covered by a torsion pendulum measuremeptpj and a 
mHz experiment at the glass transition itselieSl. Fig. 8 
shows the barrier distribution function fitting these data 
over the whole range. Again, it turns out to be possible 
to fit the a-relaxation in terms of a gaussian for the pri- 
mary barrier distribution function faiV), ascribing the 
slow rise of the damping towards the glass transition to 
secondary relaxation (the shaded area in Fig. 8). 

Fig. 9 shows the fit of the glass transition data0, cal- 
culated with Tg = 1460 K, G = 31 GPa, = 0.55 and 
Aa = 0.263 eV. If one takes the value Tg - To ^ 850 K 
and its rather large error bars from the shift factors de- 
termined inCJ, one calculates Ofg = 6 ± 4 from eq. (|3|). 
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FIG. 10. The barrier distribution function of Na20:2Si02, 
together with torsion pendulum data. References see text. 
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FIG. 11. Dynamical mechanical data of Na20:2Si02 at the 
glass transition, together with the model fit (see text). 

It is interesting to compare vitreous silica to a sodium 
silicate glass, Na20:2Si02, which has a much lower glass 
temperature Tg = 733 K and is much more fragil(£j, with 
Tg — To — 267 K. Fig. 10 shows the barrier disttibution 
function, together with torsion pendulum datac2l at Out 
Hz. The a-peak region is again fitted to the mHz dataeJ, 
shown in Fig. 11, with G — IQ-iGPa, Wa = 0.39 and 
A„ = 0.092 eV. The shift factorsEJ require = 13 ±1.5. 



IV. RESULTS AND DISCUSSION 

A. Results for the flow process 

The preceding section showed for all three examples a 
good fit of the a-relaxation in terms of a narrow gaus- 
sian primary barrier distribution function, together with 
a bigger or smaller contribution of the secondary relax- 
ation. The contribution of the secondary relaxation to 
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the a-process was biggest in the strong glass former sil- 
ica and smallest in the very fragile case of polystyrene. In 
all three cases, excellent fits of the time or frequency de- 
pendence of the a-process could be obtained (see Figs. 7, 
9 and 11). Simultaneously, the temperature dependence 
of the shift factors for could be described accurately 
by the temperature slope as of the rise of the secondary 
relaxation above Tg in eq. (|3l|). 

Table I compiles the fit parameters for the a-process 
determined in the preceding section. 



Table I: Model parameters for the a-process 



substance 


polystyrene 


Si02 


Na20:2Si02 


Tg (K) 


372 


1460 


733 


G(r,) (GPa) 


1.69 


31 


19 




0.54 


0.54 


0.39 


A„ (eV) 


0.107 


0.263 


0.092 




6.9±0.1 


6±4 


13±1.5 



In all three cases, the relative weight Wa of the a- 
process is not very far from 1/2. The a-process seems to 
occur when the secondary relaxations reduce the shear 
modulus to about half its infinite frequency value, i.e. 
when the secondary relaxational shear response equals 
the elastic one. The situation reminds one of the defi- 
nition of the Maxwell time tm = tIo/G, the time when 
the flow response equals the elastic one. Thus one could 
think of a generalized Maxwell criterion for the onset 
of the flow process, namely when the shear response by 
relaxation is the same as the elastic shear response. Intu- 
itively, this is plausible: the time at which a macroscopic 
shear stress for a given strain decays to half its initial 
value should correspond to the lifetime of the microscopic 
shear pattern of the glass former, because the microscopic 
stresses will be expected to decay on the same time scale 
as the macroscopic ones. But a decay of the microscopic 
stresses necessarily implies a coupling between different 
relaxation centers in the glass. Thus one would have to 
go beyond the simple idea of uneaupled relaxation centers 
to understand the flow processeS. 

In all three cases, the width is relatively narrow, 
less than one tenth of the barrier Va{Tg) itself, in the 
third case even less than the twentieth part. 

Finally, the coefficient as for the temperature rise of 
the secondary barrier distribution function /s(l^) above 
Tg, related to the fragility of the glass former, is not so 
very different for the three cases. There is a factor of two 
between silica and the silicate glass, but taking the dif- 
ferences in the glass temperature into account, one sees 
that the relative rise per Kelvin is the same for the two 
glasses, admittedly within large error bars. The differ- 
ence is rather in fsiV) itself; obviously, relaxing entities 
form much more easily and with much lower barriers in 
the multiply broken network of the sodium silicate glass 
than in the continuous random network of silica, making 
the glass temperature a factor of two lower. 



The large difference in fragility between silica and 
polystyrene is not so much due to a difference in as, but 
rather in the product V fs{V) at Tg (see Figures 6 and 8), 
which is about 0.3 in polystyrene and about 2.5 in silica. 
This product enters into eq. (|3|) for the Vogel-Fulcher 
temperature Tq. In this view, a substance is fragile when 
it has a low density of secondary relaxing entities at the 
relaxation time of the flow process. Naturally, the fast 
rise of the number of secondary relaxing entities above 
Tg remains a necessary condition. 

The coefficient as — 6.9 ± 0.1 for polystyrene in Table 
I, obtained from the temperature dependence of the shift 
factors (see Fig. 7), agrees within experimental error 
with the values as = 7.7 ± 1 and as = 10 ± 3 determined 
from neutron and Brillouin data above Tg (see Fig. 5). 



B. Discussion 

Here is the proper place to remind the reader that with 
all the presented equations and fits, the glass transition 
is not really explained. The Gilroy-Phillips formalism 
supplies only a reformulation of well-known equations in 
terms of a barrier distribution function. It is a new way 
to look at the data, an encouragement to seek explana- 
tions in a new direction, complementary to attemptSj-to 
understand the glass transition from the liquid sideC3, 
because it starts from a description of the relaxation in 
the glass phase. But the present work does neither ex- 
plain the fast rise of the secondary relaxation above Tg, 
nor does it supply a quantitative explanation why the 
shear modulus breaks down completely when it is halved 
by the secondary relaxation. It merely helps to quantify 
and to visualize these experimental facts. 

The concept of the barrier distribution function is 
based on the idea of independent thermally activated re- 
laxation processes in disordered surroundings. As long as 
this idea applies without any restriction, one should find 
a temperature-independent barrier distribution function, 
as one does indeed in the glass phase of polystyrene and 
silica. The quantitative comparison of measurements at 
different temperatures and frequencies with the equations 
of Section II. C enables a much more stringent check of 
the Arrhenius concept than a mere Arrhenius tempera.- 
ture shift of a broad relaxation peak with frequencyEI^EI. 
From the few such stringent checks reported so far, it 
is already clear that the concept does not always work 
perfectly well in the glass phase. In BPA-PCL-another 
amorphous polymer, there is Raman evidencal3 for an 
increase of the secondary barrier distribution function 
fs{y) with temperature below Tg , even though the 7- 
relaxation peak shift is perfeatly Arrhenius-likecfl. There 
seem to be more such casealj. This is yet another ques- 
tion calling for a closer investigation. 

The outcome of the fits of the flow process in terms 
of the primary barrier distribution function faiY) shows 
once more the collectivity of the flow process. These are 
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clearly not independent thermally activated relaxation 
processes; otherwise Va would not, shift with tempera- 
ture. There has been an attempto to explain the tem- 
perature dependence of Va (T) in terms of a proportion- 
ality to G(T) ("flow by shoving"). However, as pointed 
out in Section III.C, one finds a Griineisen Tg = 3.8 for 
G(T) in polystyrene. This is a factor of seven too weak 
to explain the temperature shift of Va {T) . 

In any case, the description of the flow process in 
terms of fa{V) in principle does not prejudice anything. 
One can absays return from faiV) to the conventional 
descriptiontj in terms of H{t) via eq. (|Tl]). What might 
be questionable is the specific subdivision into primary 
and secondary relaxation. Here, it was decided to fit the 
primary barrier distribution function fa{V) with a gaus- 
sian. On the basis of the data, this choice can be only 
justified for polystyrene; in the other two cases, one could 
choose differently. The conclusion that the shear modu- 
lus is reduced to about half its infinite frequency value 
by the secondary relaxation alone would then no longer 
be valid for these two cases. 

Independent of the choice of a specific function for 
the primary barrier distribution function, the new de- 
scription is much more convenient than the old one if 
one wants to compare different temperatures. It does 
not only allow for a stringent check of the Arrhenius be- 
haviour in the glass phase, but it also holds the promise 
to provide -a deeper understanding of the glass transi- 
tion riddleaj. As shown in Section II. E, the stretching of 
the flow process is quantitatively related to the width of 
the a-peak. More important, the Vogel-Fulcher behavior 
is an inherent feature of the generic case of a constant 
secondary barrier distribution function, provided it has 
a constant temperature rise. But one can as well un- 
derstand deviationscUa from the Vogel-Fulcher behavior 
in terms of maxima and minima of the secondary bar- 
rier distribution function. In particular, one must expect 
deviations from a single Vogel-Fulcher law when Va(T) 
sweeps throusk a secondary relaxation peak, as one in- 
deed observescH 



V. SUMMARY AND CONCLUSIONS 

The Gilroy-PhiUips modeS for secondary relaxation 
in glasses, based on the Arrhenius-Krameraj picture 
of thermally activated jumps over energy barriers, has 
been extended to describe the flow process (primary or 
a-process). The proper treatment of uncoupled ther- 
mally activated relaxation events in disordered surround- 
ings yields a temperature-independent barrier distribu- 
tion function /(V). The temperature independence holds 
as long as one has a constant number of independent re- 
laxation centers. The function has to integrate to 1 in 
order to bring the shear modulus down to zero, a very 
convenient sum rule. As it turns out, one can use the 
formalism to describe any kind of relaxation, admitting a 



temperature dependence of the barrier distribution func- 
tion. Thus one can separate the trivial Arrhenius temper- 
ature dependence from nontrivial temperature changes of 
the relaxation. 

The relaxation below the glass temperature Tg, in the 
frozen energy landscape of the glass, should be describ- 
able in terms of a temperature-independent barrier dis- 
tribution|i|mction (in some cases like poly styr ene, CKN 
and silicall3't§ it works, in others it does notllZica), but the 
flow process is most certainly of a different kind. There- 
fore one has to postulate a secondary barrier distribution 
function fs{V) for the secondary relaxation, and an ad- 
ditional primary barrier distribution function fa{V) for 
the flow process. This additional primary function shifts 
its maximum to lower values with increasing tempera- 
ture and serves simultaneously as a relatively sharp cut- 
off for the secondary relaxation. This property supplies 
a quantitative basis for the puzzling relation between the 
fragility and the fast rise of the picosecond process above 
TgL9. It explains the unusual temperature dependence of 
the flow process in terms of the temperature rise of the 
secondary relaxation, and might even shed some light on 
the detailed temperature-dependence of the flow process 
in specific substancescIrEJ. |_. 

Reformulating the classical rheological equationa£3, ex- 
perimental literature data in polystyrene, vitreous silica 
and a sodium silicate glass could be fitted in terms of the 
barrier distribution fusptions. The results corroborate 
the earlier conclusionllH of a temperature-independent 
secondary barrier distribution function in the glass phase 
of polystyrene, but show a strong increase of the sec- 
ondary relaxation with increasing temperature above Tg . 
Taking the flow process as a gaussian barrier distribu- 
tion function, the fits show the onset of the flow process 
more or less at the point where the secondary relaxation 
reduces the shear modulus to half its infinite frequency 
value. This was not only found in polystyrene, but in the 
two other cases as well. It remains to be seen whether 
such a generalized Maxwell criterion for the onset of the 
flow process is a general property of glass formers. In 
polystyrene, the temperature coefficient of the fast rise 
of fs{V) above Tg calculated from the fragility coincided 
within experimental error with the fast rise seen in neu- 
tron and Brillouin experiments. 

The findings suggest a rapid change of the energy land- 
scape above Tg, the landscape getting more and more 
rugged with more and more minima and saddle points as 
the temperature rises, thus also explaining the entropy 
crisis at the Kauzmann temperature. This does not really 
answer the central questions around the glass transition, 
but focusses the attention on two points: (i) what is the 
possible mechanism of formation of an increasing number 
of relaxing entities with increasing temperature above Tg 
(ii) why does the flow process set in when the integrated 
relaxational response to an external shear stress equals 
the elastic one. 
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Figure Captions 

Fig. 1: Asymmetric double well potential with barrier 
height V and asymmetry A as a function of a generalized 
coordinate. 

Fig. 2: Secondary barrier distribution function fs{V) 
(the shaded area), together with the cutoff by the a- 
relaxation peak (schematic). 

Fig. 3: Secondary barrier distribution function fs{V) 
calculated from literature data of the mechanical damp- 
ing of amorphous polystyrene in the glass phase at dif- 
ferent frequencies. References see text. The line is a fit; 
the same fit is also shown in Figs. 4, 5 and 6. 

Fig. 4: Secondary barrier distribution function /s(V) 
calculated from literature data of neutron and Ra- 
man scattering from amorphous polystyrene in the glass 
phase. References see text. 

Fig. 5: Secondary barrier distribution function fs{V) 
calculated from literature data of neutron and Brillouin 
data of amorphous polystyrene, both in the glass phase 
and above the glass temperature Tg. References see text. 



Fig. 6: Secondary barrier distribution function fs{V) 
calculated from literature mechanical damping and neu- 
tron data of amorphous polystyrene up to the glass tran- 
sition. The peak at the end shows the gaussian fa{V) de- 
scribing the a- process at the glass temperature Tg. Ref- 
erences see text. 

Fig. 7: Creep data of the a-process in polystyrene at 
the glass transition (reference see text), together with a 
fit in terms of f{V). 

Fig. 8: The barrier distribution function of vitreous sil- 
ica, together with torsion pendulum and Brillouin damp- 
ing data. References see text. 

Fig. 9: Dynamical mechanical data of vitreous silica 
at the glass transition, together with the model fit (see 
text). 

Fig. 10: The barrier distribution function of 
Na20:2Si02, together with torsion pendulum data. Ref- 
erences see text. 

Fig. 11: Dynamical mechanical data of Na2 0:2Si02 
at the glass transition, together with the model fit (see 
text). 
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